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We investigate Hamiltonian systems with two degrees of freedom by 
using renormalization group method. We show that the original Hamil- 
i 1 tonian systems and the renormalization group equations are integrable if 

the renormalization group equations are Hamiltonian systems up to the 
c~| | second leading order of a small parameter. 

O 

Many systems in the nature are described by Hamiltonian equations of motion, 
and we are interested in dynamical properties of the systems. Although behavior 
of systems is revealed by solving equations of motion, we cannot generally obtain 
exact solutions to Hamiltonian systems with more than one degree of freedom 
because of lack of integrals. We therefore solve equations of motion by using 
perturbation techniques. Naive perturbation often fails to give global solutions 
since it yields secular terms, and hence improved perturbation methods have 
been developed, for instance, averaging method, multiple scale method, matched 
asymptotic expansions and canonical perturbation theory [^]. 

Recently, renormalization group method [| is proposed as one of the most 
powerful tools to construct approximate solutions since it unifies many of the 
perturbation techniques listed above. This method reduces equations of mo- 
tion to amplitude equations called renormalization group equations (RGEs) by 
ignoring fast motion. The reduced equations, RGEs, must reflect features of orig- 
inal systems, and we expect that dynamical properties of the original systems are 
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obtained from RGEs. In this article, we consider two types of perturbed Hamilto- 
nian systems with two degrees of freedom. We investigate symplectic properties, 
integrability and integrals of RGEs, which are related to properties in the original 
systems. We first present our main result: 

Theorem Let Hamiltonian be represented as follows: 

H(qi,q 2 ,Pi,p 2 ) = Ho(q l ,q 2 ,Pi,p 2 ) + eV l (qi,q 2 ), (1) 

#0(91,92, pi, P2) = + + qf + q 2 ), ( 2 ) 

where the potential Vi(qi, q 2 ) is a homogeneous cubic or quartic function of q\, q 2 
and e is a small parameter, i.e. |e| << 1. If renormalization group equation of 
the system ([![) is a Hamiltonian system up to the second leading order of e, then 
the original system ^) and the renormalization group equation are integrable. 

Let us briefly review renormalization group method by using a simple system 
with one degree of freedom: 

H(q,p) = H {q,p)+eV 1 (q), (3) 
H (q,p) = \(p 2 + q 2 ), V l {q) = \q 2 - (4) 



The equation of motion is 



W +q = ^ (5) 



and the exact solution is 

q = B cos(v / l + et) + C sin(v / l + et), (6) 

where Bo and Co are constants of integration and determined by initial condition 
at the initial time t = t . We perturbatively solve Eq.@ by expanding q as a 
series of positive powers of e: 

q = q + eq 1 + e 2 q 2 H . (7) 

This naive expansion gives 

q(t; t , -Bo, Co) = B Q cost + Co sint + - (t — to) (Co cost — B sint) 



e 2 



(t- t ) (C cost - B sint) + (t - t ) 2 (B cost + C sint) +0(e 3 ), (8 



and this expansion breaks when e(t — to) > 1 because of secular terms. Here we 
ignored homogeneous parts of qj(j > 1) which is the kernel of the linear operator 
L = d 2 /dt 2 + 1, since they can be included in q . Renormalization group method 
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removes the secular terms by regarding Bq and Co as functions of the initial time 
t , and the evolutions of B(t ) and C(t ) are determined by RGE || f3|: 



dq 



= 0, for any t. (9) 

t =t 



In other words, the RGE is 

)C(t)+0(e 3 ), (10a) 
5(t) + 0(e 3 ), (10b) 



dB(t) l e r 

dt \2 
dC(t) ft e 2 " 



dt \2 ; 

and it gives the renormalized solution q 
q RG (t) = q(t;t,B(t),C(t)) 



e e 2 \ \ „ // e e 2> 



= S cos 1(1 + - - — I tj + C sin Ml + - - — I tj + 0(e 3 ),(ll) 

where B = B(to) and Co = C(to). This expression gives an approximate but 
global solution to Eq.(||) up to 0(e 2 ). 

In the RGE (|1(]) B(t) and C(t) are canonical conjugate variables because the 
equation is yielded by Hamiltonian 



H RG =( € --^)(B 2 + C 2 ), (12) 



and 

dB dH RG 



dt dC 



+ 0(e 3 ), (13a) 



Moreover, the system governed by the H RG is obviously integrable and (B 2 + 
C 2 )/2, which corresponds to H , is an integral of H RG . Do these statements hold 
even in systems where chaotic motion appears? 

To answer this question, we calculate RGEs in system ([]]), whose equation of 
motion is 

First let potential V\ be homogeneous quartic functions of qi,q2- 

Vi(qi, Qz) = ot\q\ + cn 2 qfq2 + a 3 qfqj + a 4 <?i<?2 + a. h q\. (15) 
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We expand q\ and q 2 as 

Qj = Qjfi + e ?i,i + e2( lj,2 + •'•> 0' = 1> 2), (16) 
and write the zero-th order solution to the equation (|14]) as 

q ji0 = B 3 cos(t) + Cj sin(t), (j = 1, 2), (17) 

where we omit subscript of Bj and Cj (j = 1,2) for simplicity of symbols. 
Following the procedure to obtain RGE, we get 

dB 

= fj( B u Ci, B 2 , C 2 ; e, a u • • • , a 5 ), (18a) 

dC 

-^f = 9j( B i, Ci, B 2 , C 2 ; e, a 1} ■ ■ ■ , a 5 ), (18b) 

although the explicit forms of fj and gj are not shown here because of complexity. 

To clarify the condition of «i, • • • , with which the RGEs flT8j ) become Hamil- 
tonian systems, we define Aj as 

For a RGE, both of the equations Ai = and A 2 = are satisfied if and only if 
Hamiltonian exists and 

dB j dH RG . . 

(20a) 



(20b) 



dt dCj ' 

dCj dH RG 
~dT ~ dBj 

We show the Aj up to the second order of e: 

,2 

Ai = -(B 1 C 2 -C 1 B 2 ){(9a 2 2 + 4al-24a 1 a 3 -9a 2 a 4 )(B 1 B 2 + C 1 C 2 ) 
8 

+ 3[(a 2 + a 4 )a 3 - 6(ai« 4 + a 2 a 5 )](B% + C|)| , (21a) 

2 

A 2 = --(5 1 C , 2 -Ci5 2 ){(9a2 + 4a2-24a3a 5 -9a 2 a 4 )( J B 1J B 2 + CiC 2 ) 
8 

+ 3[(a 2 + a 4 )a 3 -6(aia 4 + a 2 a 5 )](E 2 + C' 2 )}. (21b) 

The Ai and A 2 take zeros irrespective of values of Bj and Cj if and only if the 
three equations hold: 

9a 2 , + 4a 2 - 24ai« 3 - 9a 2 a 4 = 0, (22a) 
9a 2 + 4a 2 - 24a 3 a 5 - 9a 2 a 4 = 0, (22b) 
(a 2 + a 4 )a 3 — 6(ai« 4 + a 2 «s) = 0. (22c) 
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Consequently, RGEs are Hamiltonian systems when the coefficients of V x , namely 
«i, • • • , «5, satisfy the condition ( p2|) . We remark that Aj has no terms of order 
e, which is the leading order of Eq. (p"8|) , and hence the RGEs are Hamiltonian 
systems for any coefficients up to 0(e). 

To understand the meaning of the condition (E^), we use the Bertrand- 
Darboux theorem||. with which we judge integrability of "natural" Hamiltonian 
systems with two degrees of freedom whose form is 

H=±(£+pt) + V(x,y). (23) 

The theorem states that the following three conditions are equivalent: 

(1) There is an integral which is independent of Hamiltonian and quadratic 
with respect to the momenta. 

(2) For a set of constants, (a, b, b', c±, c 2 ) ^ (0, 0, 0, 0, 0), the potential V(x, y) 
satisfies the Darboux equation 

( v m ~ v xx)(-2axy ~ b'y - bx + c x ) 
+2V xy (ay 2 — ax 2 + by — b'x + c 2 ) 

+V x (6ay + 3b) + V y (-6ax - 3b') = 0, (24) 

where V x = dV/dx, V xx = d 2 V/dx 2 and so on. 

(3) The system is separable in Cartesian, polar, elliptic or parabolic coordi- 
nates. 



We apply the Bertrand-Darboux theorem to the original Hamiltonian sys- 
tems written by Eqs.(|l|) and flT5|) to judge whether the systems are integrable 
or not. The Darboux equation restricts the coefficients a±, ■ ■ ■ , as in the poten- 
tial function ( |T5D (Table I), and we can find in the Table [I] the condition fl2"2|). 
Consequently, when RGEs become Hamiltonian systems, the original systems 
are integrable. Moreover, the RGEs are also integrable since they have another 
integral J, 

/ = (Bf + C\ + Bi + CD/2, (25) 

except for Hamiltonian H RG itself. Note the quantity / corresponds to Ho in 
Eq-®- 

For the potential of homogeneous cubic functions of q x and q 2 , 

Vi(q) = a x ql + a 2 qfq 2 + a 3 <?i<?2 + "4^2, (26) 

secular terms appear only in even orders of e and the leading and the second 
leading orders of RGEs are e 2 and e 4 , respectively. Then we calculate A x and A 2 
up to the fourth order of e and 

A x = —{BiC 2 - B 2 Ci)(3(aia 3 + a 2 a±) - a 2 - (27a) 
A 2 = -w( B i C, 2-5 2 C 1 )(3(a 1 a 3 + a 2 «4) -a 2 2 -aj)h 2 , (27b) 
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Table 1: Condition of ax, ■ ■ ■ , as, to satisfy the Darboux equation. Here ax, • • ■ , as 
are the coefficients of homogeneous quartic potential functions. There are two 
cases I and II, and case II must simultaneously satisfy the three equations, which 
coincide with the condition that RGEs become Hamiltonian systems. 





Condition for aj 


Potential function 


I 


ot2 = Q-4 = 0, 03 = 2ax = 2a§ 


V 1 = a l {ql + qlf 


II 


9o| + 4a| — 24oio 3 — 9o 2 o 4 = 
9o| + 4«3 - 24a 3 a 5 - 9o 2 o 4 = 
(a 2 + 0:4)03 — 6(0404 + 02O5) = 


Vx = axqf + a 2 q\q 2 + a 3 q\q\ 
+a 4 qxq 2 + a^q\ 



where 

hx = 40a 2 (3a 1 + a 3 )(B 2 1 +Cl) 

+ (66oia 2 + 139o 2 a 3 — 495oi0 4 + 1860304) [B\ + C 2 ) 

+ 6(33o? + 34o^ + 12o^-40oio 3 + 15o2O4)( J Bi J B 2 + CiC2), (28a) 

h 2 = 40o 3 (3o 4 + o 2 )(5 2 2 + C , 2 2 ) 

+ (6603O4 + I3902O3 — 495oi04 + 186oio 2 )(i? 2 + Cl) 

+ 6(33o 2 + 34o 2 + 12o 2 -40o 2 O4 + 15oiO3)( J Bi J B2 + ^iC 2 ). (28b) 

One condition for Ai = A 2 = is hx = h 2 = whose unique solution is the 
trivial one, i.e. (01, o 2 , o 3 , o 4 ) = (0,0,0,0), and hence we ignore this condition. 
The other condition is 3(oi0 3 + o 2 o 4 ) — a\ — a\ = 0, and this condition is found 
in Table which shows the condition for the Darboux equation being satisfied. 
Consequently, the original systems (P with the cubic perturbation are integrable 
when RGEs are Hamiltonian systems up to the second leading order. We therefore 
proved our Theorem. 

In summary, we investigated dynamical properties of renormalization group 
equation (RGE) which are symplectic properties, integrability and integrals. 
Hamiltonian systems with two degrees of freedom were considered, whose inte- 
grable part is harmonic oscillators and homogeneous cubic and quartic potential 
functions are added as perturbation. The main result of this article is that the 
original Hamiltonian systems and RGEs are integrable if the RGEs are Hamilto- 
nian systems up to the second leading order of a small parameter. 

We have two future works. One is generalization of potential function. We 
can directly compare the Darboux equation with condition for RGE being Hamil- 
tonian systems if we write RGEs for the general potential. The other is extension 
to systems with higher degrees of freedom. For the purpose, extended Bertrand- 
Darboux theorem is available. Taking the contraposition of our main result, 
we suppose that we can obtain information on chaotic behavior by considering 
how symplectic properties of RGEs break in non-integrable systems. 
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Table 2: Condition of cci, • • • , 0:4, to satisfy the Darboux equation. Here oti, ■ ■ ■ , a± 
are the coefficients of homogeneous cubic potential functions. Four cases I,- • •, IV 
are obtained, and the condition in case II coincides to the one that RGEs become 
Hamiltonian systems. 





Condition of otj 


Potential function 


I 

II 
III 

IV 


Oil — OL2 = «3 = OiA = 

3(0:10:3 + 0:2*^4) — ot\ — a\ = 

OL\ = 2o;3, (X2 = 0:4 = 

cti = «3 = 0, 2a 2 = 0:4 


Vi = 
Vi = otxql + a 2 qfq2 

+"3^1 ql + oi^q\ 
V x = a 3 (2qf + qi ql) 
V x =a 2 {2q\ + qfq 2 ) 
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